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S ■ ABSTRACT 
a ; 

^ ^ , Obliquity tides are a potentially important source of heat for extrasolar planets on 

Q I close-in orbits. Although tidal dissipation will usually reduce the obliquity to zero, a 

CSJ I nonzero obliquity can persist if the planet is in a Cassini state, a resonance between spin 

^ ' precession and orbital precession. Obliquity tides might be the cause of the anomalously 

^> ■ large size of the transiting planet HD 209458b. 

oo ; 
\o . 

Tij- . Subject headings: stars: individual (HD 209458) — planets and satellites — planetary sys- 

. terns: formation 

in : 

O ; 

' 1. Introduction 

2 . A wonderful surprise of the last decade was the discovery of "hot Jupiters" : giant planets that 

'• orbit their parent stars at a distance smaller than 10 stellar radii (Mayor & Queloz 1995, Butler et 

I al. 1997). The proximity of the star and planet raises the possibility of tidal interactions between 

. ^ I them. Of these interactions, the most frequently discussed are the synchronization of the planetary 

^ ' rotation period and orbital period, which happens over ~10^ yr, and orbital circularization, over 

. 10®^^ yr (see, e.g., Rasio et al. 1996, Lin & Gu 2004). Less frequently discussed, but often implicit, 

is tidal evolution of the obliquity (the angle between the planetary spin axis and the orbit normal), 
which occurs on the same time scale as synchronization (Peale 1999). 

We wish to point out that in some circumstances a close-in planet can maintain a nonzero 
obliquity, and that the consequent heat from obliquity tides can be large enough to affect the 
internal structure of the planet. We present the heating calculation first, in § 1, since it does not 
depend upon our particular scheme for maintaining the obliquity. Our proposal, given in § 2, is 
that hot Jupiters occupy Cassini states, in which spin precession resonates with orbital precession. 
In § 3 we ask whether obliquity tides can explain the famously small density of the extrasolar planet 
HD 209458b. In § 4, we discuss the implications, strengths, weaknesses, and possible tests of this 
theory. 
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2. The power of obliquity tides 



Tidal torques synchronize the planetary spin frequency (lu) and orbital mean motion (n), and 
often reduce the obliquity (9) to zero (Goldreich & Peale 1970). Henceforth, if the orbit is circular, 
the tidal bulge is motionless in the reference frame of the planet and energy dissipation ceases. 
However, dissipation continues if the planet somehow maintains a nonzero eccentricity or obliquity. 
For a star of mass My, and a planet of mass Mp and radius Rp in an orbit with semimajor axis a, 
the rate of energy loss through tidal friction is 



dE 



9 h fGMl 

Tl I 

log V 



(7e2 + sin^ 9) . 



(1) 



This formula, derived in a different context by Wisdom (2004), employs the customary model 
of friction within astronomical bodies: Q is the "quality factor" of tidal oscillations (the inverse 
of the fractional energy dissipated per cycle). The factor h is the displacement Love number, 
parameterizing our ignorance of the planet's deformability. The factor in square brackets may be 
written 

106 V 3 days J Vi?jup 



2 X 10^^ erg s"^ 
where P is the orbital period. 



(2) 



For this heat source to play a significant role in determining the structure of the planet, it 
must be comparable to the intrinsic luminosity Lq of other processes (e.g. gravitational contraction, 
stellar irradiation, radioactivity, and nuclear reactions). For Jupiter, Lq = 3 x 10^^ erg For 
hot Jupiters, theoretical estimates of Lq are difficult to summarize (depending as they do upon the 
planet's age, temperature, composition, atmosphere, and interior structure) but for billion-year- 
old planets of Jupiter's radius and mass, one expects Lq ~ 2 x 10^"'' erg (see, e.g., Guillot & 
Showman 2002, Baraffe et al. 2003). Tidal heating is dominant when e ^ 0.04 oi 9 ^ 0.1. The 
importance of eccentricity tides is well known (Bodenheimer, Lin, k. Mardling 2001), but the 
equivalent importance of obliquity tides does not seem to have been appreciated. 



3. Cassini states for extrasolar planets 

The reason is probably a mistaken intuition that = is the only possible endpoint of tidal 
evolution. In general, the planet's spin axis and its orbit both precess, in response to additional 
planets, satellites, the stellar quadrupole, or other torquing agents. With spin and orbital preces- 
sion, the outcomes of tidal evolution are Cassini states (Colombo 1966, Peale 1969, Ward 1975), 
in which the orbit normal fi and spin axis s precess at the same rate about the same axis k (see 
Fig. 1). For a given body, there are at most two stable Cassini states, differing in whether s and 
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n are on the same side of k (state 1) or opposite sides (state 2).-'^ The obhquity 9 = cos~^(s • n) 

of a Cassini state is not generally zero, nor is it necessarily close to zero. Although many satellites 
in the Solar system are in state 1 with vanishingly small obliquity, the Moon is in state 2 with 
e = 6?5.2 




Planet 



Star 




Fig. 1. — (Left.) Illustration of the spin axis s, orbit normal n, precession axis k, orbital inclination 
/ and obliquity 6. A general orientation for s is depicted. In a Cassini state, s is coplanar with k 
and n. (Right.) A convenient set of Cartesian axes, in which n is along the Z axis, k is in the YZ 
plane, and (^,'0) are traditional polar angles describing s. In a Cassini state, = ±90°. 

Idealizing the planet as an oblate rigid body on a circular and uniformly precessing orbit of 
inclination I = cos~-^(k • n), the Cassini obliquities 0i obey 

cos 6i sin 9i — e s\n{6i — /) = 0. (3) 

Here we have defined e = —g/a^ where g is the nodal precession frequency and a is the spin 
precessional constant for a fixed orbit. The latter can be written 
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where C > B = A are the planet's principal moments of inertia, and u is its spin frequency. 
During synchronization, a; — n. When e < Ccrit = 

(sin2/^7 + cos2/3 7)-3/2^ Eq. (3) has 4 roots, 
corresponding to the two stable states, 1 and 2, and two unstable states, 3 and 4. For e > ecrit) 
there are only 2 roots, corresponding to states 2 and 3. Eq. (3) can be derived from the governing 
Hamiltonian under the assumption of principal-axis rotation (see, e.g.. Ward 1975), 

a 



n = no-^{h-sy-g{k-s), (5) 

where Hq is the s-independent portion. Fig. 2 shows contours of H{s) for an illustrative case. 



^Formally, there are two other equihbria: state 3 is hneaxly stable but is unstable to tidal evolution (Goldreich & 
Peale 1970), and state 4 is linearly unstable. 

^The observation that the Moon's s, fi, and k are coplanar was the basis of G. D. Cassini's third law of lunar 
motion, formulated in 1693. 
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Fig. 2. — (Left.) Contours of 7i(s), for e = 0.6 and / = 0.1. The unit sphere is projected onto the 
XY plane. The thick contour is the separatrix, which divides the basins of attraction of states 1 
and 2 under tidal evolution. The numbered dots show states 1, 2, and 4. The retrograde state 3 
is not shown. The separatrix angles 9s are marked by squares. (Right.) Cassini state obliquities 6i 
(solid lines) and separatrix angles 6s (dashed lines) as a function of e, for 7 = 0.1. 

A hot Jupiter settles into a Cassini state on the same ~10^ yr time scale as spin-orbit syn- 
chronization. Whether it ends up in state 1 or 2 depends upon its initial obliquity and orbital 
inclination, and upon e, which may be a function of time. Reasons for variations in e include the 
disappearance of the protoplanetary disk, migration of the planet or its fellow planets or satellites, 
contraction of the planet, and spin alteration of the planet or star. If e varies slowly compared to 
g''^ and ( "adiabatically" ) , and slowly compared to the Cassini state settling time, then the 
obliquity tracks the evolving Cassini obliquity. Abrupt changes in e or may cause the planet to 
leave a Cassini state, and ultimately to switch states if it lands in the basin of attraction of the 
other stable state. 

Fig. 2 suggests three ways in which a hot Jupiter can maintain a significant obliquity. First, if 
e > Ecrit, the only stable state is 2, for which the obliquity is nonzero (02 — as e — oo). This is 
the case for the Moon. Second, it could be in state 1 or 2 with e ~ ecrit, but this requires a special 
coincidence, since e involves both planet-specific properties (namely, its moments of inertia) and 
seemingly unrelated quantities specific to other bodies (which determine the precessional torques). 
Third, the planet could be in state 2 with e < ecrit- When e is small, the basin of attraction of 
state 2 is also small, but it is possible to be captured into state 2 when e is large and then tipped 
over to large obliquity as e is gradually reduced. Ward & Hamilton (2004) proposed this scenario 
to explain why Saturn has a much larger obliquity than Jupiter. In the next section, we propose a 
similar scenario for a particular extrasolar planet. 
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4. The case of HD 209458b 

The planet HD 209458b transits its parent star (Charbonneau et al. 2000, Henry et al. 2000), a 
fortuitous circumstance that enables many interesting measurements, including that of the planet's 
mean density, which is 0.33 g cm~^. This is 27% of the Jovian value and is the smallest mean 
density of all 7 known transiting extrasolar planets (see, e.g., Alonso ct al. 2004, Konacki et al. 
2005, Pont et al. 2005). Theorists have struggled to explain this anomaly, usually by attempting to 
identify an overlooked internal heat source, although Burrows, Sudarsky, & Hubbard (2003) argued 
that the density is not so terribly anomalous. We present a new hypothesis: HD 209458b resides 
in Cassini state 2 with a large obliquity, whereas most hot Jupiters reside in state 1 with small 
obliquities. 

We are led to imagine the following sequence of events: (1) The planet forms at a large orbital 
distance, with a nonzero (but not necessarily large) obliquity. (2) The planet migrates inward to 
its current position. (3) As the spin and orbit are synchronized over ^^10^ yr, the planet falls 
into Cassini state 2, whether by chance or because state 2 is the only possibility. (4) As the disk 
disappears over ~10^ yr, the orbital precession rate (and e) decreases, forcing the obliquity to grow. 
(5) The planet remains in state 2 for billions of years. 

Is this scenario compatible with order-of-magnitude estimates of the relevant quantities? We 
assume I = 0.1, because Winn et al. (2005) found an angle of «4° between the sky projections 
of the stellar spin axis and orbit normal, and it seems reasonable that / is of the same order of 
magnitude. We estimate (C — A)/C by scaling Jupiter's observed value (de Pater & Lissauer 
2001) by (a;/a;jup)^. Application of Eq. (4) gives a/n ^ 10~^ after synchronization. We assume 
that g was initially determined by the torque of the protoplanetary disk, which we idealize as a 
minimum-mass solar nebula. The precession rate is sensitive to the (unknown) inner radius of the 
disk: for rmm = 1.6a, \g/n\ = 5 x 10~^, rising to 10""^ for 1.3a. These estimates ignore resonance 
effects, which we find through trial numerical integrations to be capable of enhancing g by an order 
of magnitude. As the disk disappears, g falls until other torques dominate, such as those from 
the stellar quadrupole or additional orbiting bodies. In the former case, assuming J2 = 10~^, the 
precession slows to \g/n\ = 4 x 10"^. There is no evidence for additional bodies in current radial 
velocity data (Laughlin et al. 2005), but even a body too small or too distant to have been detected 
could nevertheless cause faster precession than the stellar quadrupole alone. For example, a IM^ 
planet in a 1.2 day orbit would cause \gln\ = 2 x 10~^, while producing undetectable radial velocity 
variations of ~0.5 m s~^. We note that g is measurable, in principle, through time variations of 
the transit duration (Miralda-Escude 2002). 

With these considerations, it seems plausible that e had an initial value ranging from 0.02 
to near unity, and a final value between 10~^ and 10~^. The initial value of e could have been 
large enough that capture into Cassini state 2 occurred with reasonable probability, or even 100% 
probability if e > Ccrit- The final value of e is small enough to force ^2 90°. Thus, the current 
obliquity must be quite large. 
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The energy dissipated through obhquity tides comes at the expense of the orbital energy. We 
must examine whether it is possible for the dissipation rate to be large enough to inflate the planet, 
yet slow enough to allow the orbit to survive for 5 x 10^ yr, the approximate main-sequence age 
of the star (Cody & Sasselov 2002). To inflate the planet, Bodenheimer, Laughlin, & Lin (2003) 
found the required power to be 4 x 10^^ erg if the planet has a dense core, and an order of 
magnitude smaller if the planet is coreless. With reference to Eq. (1), this corresponds to an upper 
limit on Q/h of 10^ for a cored planet and 10^ for a coreless planet. The condition 

a GM^Mp/2a ^ ^ ^ . 



\da/dt\ dE/dt 



> 5 X 10^ yr, (6) 



yields a lower limit on of 5 x 10^. In reality, Eq. (6) is probably too restrictive by a factor of 
a few; the dissipation rate was smaller in the past, when the planet had a smaller radius, a larger 
orbital distance, and possibly a larger mass (Vidal-Madjar et al. 2003). Thus, Q/h should be ~10^ 
if the planet has a core and 10^-10^ if it does not. In comparison, for Jupiter it is thought that 
Q is between 10^ and 10® (Goldreich & Soter 1966) and h ^ 0.6 (Gavrilov & Zharkov 1977). The 
required Q/h for HD 209458b is comparable to, or somewhat larger than, the nominal Jovian value. 

The last chapter of the story is that the planet remains in Cassini state 2 for billions of years. 
Here the difficulty is that as e decreases, the width of the resonance also decreases, reducing the 
robustness of state 2 to perturbations. For e = 10~^, the angle between the separatrix angles (see 
Fig. 2) is only 0?2. For e = 10~^, it is increased to 7°. The impact of a body of mass m at escape 
velocity would produce a maximum obliquity shift of 

mvescRp my/2GMpRp ^ ^qo f 



Cn 

where the maximum is achieved for grazing incidence at the pole. Hence we must also suppose 
HD 209458b suffered no major collisions after the disappearance of the protoplanetary disk. 



5. Discussion 

Hot Jupiters should be in Cassini states. Since the obliquity of a Cassini state is not necessarily 

small, obliquity tides arc a potentially important internal heat source. Whether a given planet can 
maintain a significant obliquity depends upon its initial obliquity as well as its precessional and 
collisional histories. Obliquity tides could be the "missing" heat source that bloats the transiting 
planet HD 209458b. The implications of this hypothesis are that the planet's obliquity is nearly 
90°, its tidal dissipation factor Q and displacement Love number h obey 10® < Q/h < 10^, and 
it had a quiescent history of no major collisions after its orbital precession rate declined from an 
initially large value. 

This hypothesis has certain strengths and weaknesses relative to the two leading prior hy- 
potheses. Guillot & Showman (2002) proposed that atmospheric circulation patterns convert a 
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small fraction of the stellar radiation into heat deep within the planet. An advantage of this hy- 
pothesis over ours is that no fine tuning of Q/h or the collision history is required. A disadvantage 
is that it is not obvious why the other hot Jupiters should not experience the same phenomenon. In 
contrast, Cassini state 2 is naturally a minority state for hot Jupiters. After landing in state 2 (pos- 
sibly by chance), a planet must avoid being jostled into state 1 by collisions or the synchronization 
process. 

In the scenario of Bodenheimer et al. (2001), tidal heating from orbital circularization is ex- 
tended indefinitely because of a periodic eccentricity exchange with another planet. The extra heat 
is the fault of a well-placed third body, which is naturally expected to be a rare occurrence. Un- 
fortunately, no such body has yet been detected, and the current eccentricity is small (Laughlin et 
al. 2005, Winn et al. 2005). These problems have made the hypothesis less appealing than it once 
was, though they might be overcome through tuning of Q/h or other embellishments. A relative 
strength of our hypothesis that no third body is required, although a third body could increase the 
robustness of the high-obliquity state. 

Further work is needed to improve upon our order-of-magnitude calculations. The synchro- 
nization process must be followed in detail because a changes on the same time scale as obliquity 
evolution, and is therefore non-adiabatic (Peale 1974). We treated the planet as an oblate spheroid, 
but in reality it is triaxial due to tidal distortion, a complexity that alters the Cassini obliquities 
and dynamics (Peale 1969). We ignored any evolution of e or /, and any non-uniform precession. 
Finally, a realistic description of perturbations is needed to estimate the capture probability and 
lifetime of Cassini state 2. 

We can think of two possible tests of the theory that HD 209458b has a large obliquity, 
neither of which is easy. First, the obliquity can be measured or bounded through high-precision 
photometry of the transit (Seager &; Hui 2002, Barnes k, Fortney 2003), but the expected signal is 
smaller than 10~^ in relative flux. Second, a hot Jupiter with 9 90° has no permanent dayside or 
nightside, and should have a smaller day-night temperature difference than a planet with ^ = 0°. 
To quantify the expected temperature difference, atmospheric models of hot Jupiters (Showman 
k Guillot 2002, Cho et al. 2003, Menou et al. 2003, Burkert et al. 2005) should be generalized to 
cases of large obliquity. 

We are indebted to J. Wisdom for his suggestion to investigate Cassini states. We thank R. 
Crocker, S. Gaudi, G. Laughlin, D. Lin, S. Peale, K. Penev, D. Sasselov, K. Stanek, and W. Ward for 
helpful discussions. Work by J.N.W. was supported by NASA through grant HST-HF-01180.02-A. 
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